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In this work we consider the possibility of describing the current evolution of the universe, without 
the introduction of any cosmological constant or dark energy (DE), by modifying the Einstein-Hilbert 
(EH) action. In the context of the f(R) gravities within the metric formalism, we show that it is 
possible to find an action without cosmological constant which exactly reproduces the behavior of the 
EH action with cosmological constant. In addition the f(R) action is analytical at the origin having 
Minkowski and Schwarzschild as vacuum solutions. The found f(R) action is highly non-trivial and 
must be written in terms of hypergeometric functions but, in spite of looking somewhat artificial, 
it shows that the cosmological constant, or more generally the DE, is not a logical necessity. 
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PACS numbers: 98.80.-k, 04.50.+h 

One of the most important recent scientific discoveries 
is the accelerated expansion of the universe. Different 
data from type la supernovae 0, 0, Q observation, large 
structure information and delicate measurements of the 
cosmic microwave background (CMB) anisotropics (par- 
ticularly those from the Wilkinson Microwave Anisotropy 
Probe (WMAP) 0) have concluded that our universe is 
expanding at an increasing rate. This fact sets the very 
urgent problem of finding the cause for this speed-up. 

Usual explanations belong to one of the following three 
classes: First one reconciles this acceleration with Gen- 
eral Relativity (GR) by invoking a strange cosmic fluid, 
DE, with a state equation p = cop where u> is very close 
to -1, i.e. the fluid has a large negative pressure. For the 
particular case w = — 1 this fluid behaves just as a cos- 
mological constant A. Within this approach recent data 
obtained by WMAP correspond to the cosmological pa- 
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H (t = t ) = lOOhkms-^Ipc- 1 with h = 0.71 
to = Uoday The main problem of this kind of description 
is that the fitted A value seems to be about 55 orders of 
magnitude smaller than expected (the cosmological con- 
stant problem). The second type of explanations con- 
sider a dynamical DE by introducing a new scalar field. 
Finally the third one is trying to explain the cosmic ac- 
celeration as a consequence of new gravitational physics 
. EH action modifications have been widely considered 
in the literature 0, 0, , firstly to describe inflation, and 
more recently to describe the current cosmic speed-up, 
or even both cosmological eras simultaneously. 

The simplest way of modifying EH action is by adding 
some function f(R) with the required properties (f(R) 
gravities). For example in || it was introduced a gravita- 
tional model where f(R) = — /i 4 /i?, being the total grav- 
itational action proportional to R — /i 4 / R. This proposal 



has very interesting cosmological properties and triggered 
a lot of work on f(R) gravities applied to cosmology. 
However this kind of actions with negative powers of the 
curvature has the very serious drawback of not having 
vacuum solutions with vanishing curvature. For instance 
in the mentioned model the vacuum constant curvature 
solution is R — ±^/3/i 2 . Thus, even if one succeeds in re- 
producing cosmic acceleration, paradigmatic GR vacuum 
solutions assumed to play a major role in any fundamen- 
tal theory of gravity, such as Minkowski or Schwarzschild, 
are excluded. Other f(R) functions recently considered 
in the literature face similar problems and moreover could 
be in conflict with Solar System experiments 0, 0] . 

In this work we address the issue of finding a f(R) grav- 
ity able to reproduce the current cosmic speed-up with- 
out any cosmological constant but having R = as vac- 
uum solution. From a more formal point of view we are 
seeking for a f{R) gravity having the same Friedmann- 
Robertson- Walker (FRW) solution as the standard EH 
action with cosmological constant for non relativistic 
matter (p — 0), but being analytical at R = 0. Clearly 
the f(R) expansion at R — must start at the R 2 term 
to avoid having cosmological constant or to redefine the 
Newton constant. 

In order to consider such as the standard EH cosmo- 
logical solution with DE as f(R) gravity cosmological 
solution without DE in the same setting, we start from 
a general action S = Sg + Sm + Sde where Sq is the 
gravitational action given by: 
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d 4 x^/\g\[R + f(R)] 
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with k = 87tGjv and f(R) being any arbitrary function of 
the scalar curvature R = g^R^, where the Ricci tensor 

and the curvature tensor is de- 



where T % are 
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is given by i? M „ = R 
fined with the convention Rg^ 
the symbols for the Levi-Civita connection since we are 
using the metric formalism (see El for a recent review 
on non- metric formalisms). In addition, Sm is the ac- 
tion describing matter (also including dark matter) and 
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Sde is the action describing DE, which in particular in- 
cludes any possible cosmological constant A. For these 
two actions the corresponding energy-momentum tensor 
are given by = 2\g\- 1 / 2 SS x /Sg^ with X = M (mat- 
ter) or X = DE. For the sake of simplicity DE will be 
assumed to follow pr>E = — Pd e as state equation (i.e. 
it is just a cosmological constant) where P de = A/«. 
Thus in our notation (TdeY v = ~PdeS m v - Assuming 
that matter (including dark matter) can be described as 
a perfect fluid, the corresponding energy-momentum ten- 
sor is [T M y v = - diag(p M , ~Pm, ~Pm, ~Pm)- 

In the metric formalism field tensorial equations are 
found by performing variations of the above action 1)1)1 
with respect to the metric. Thus the equations are: 

[l+f'(R)]R^-):[R+f{R)]9^+V^f'{R) = kT^ (2) 



where ' represents the derivative with respect to R, 
P M „ = D^D U — g^ v □, □ = D a D a and D is the usu- 
al covariant derivative. By computing covariant deriva- 
tive of we find the equations of motion D^T^ = 
independently from f(R) and A. Clearly for the particu- 
lar case f(R) = we recover the standard Einstein field 
equations with cosmological constant, i.e. 



(3) 



and the matter equation of motion D^Tm)^ = 0. In 
the following we will be interested in the cosmological 
solutions of the above equations with flat spatial sections. 
Thus we will consider the line element 
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From this metric the matter equation of motion in the 
general case reads: 



pM + 3(1 + w m )pm- = 



(5) 



with dkPM = where k runs through r, 9 and <j), the 
dot represents the time derivative and we have assumed 
the matter state equation to be pm = ^mPm- These two 
equations can be integrated to give: 



PM(t) = Pm(*o) 
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where to is the present time. Now from the //, v = t 
component of we have: 

R t t [l + f'(R)]-^{R + f(R)}+A + V\f'(R) = -k Pm . 

(7) 

In a flat universe it is possible to write the scalar curva- 
ture in terms of scale parameter a = a(t) and its deriva- 
tives as: 
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R = 6 



(8) 



and subsequently, for those universes R l t = 3a/a and 
V t J'(R) = -MRf"(R)/a. The most recent cosmo- 
logical data quoted in the introduction are compatible 
with a cosmological model based on a flat FRW met- 
ric like iQJ together with Einstein field equations with 
cosmological constant A ^ and non relativistic (dust) 
matter (including dark matter), i.e. pu = 0. In this 
case we can use the equation J5J with f(R) = and 
T% = - diag(pM + PDE,-pDE,-pDE,~pDE)- Then we 
find: 
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where the subindex means that we are using stan- 
dard EH cosmology equations with cosmological con- 
stant. This notation will be relevant later on when we 
will compare standard cosmology with the results coming 
from the f(R) action for gravity that we have found in 
this work. The above equation can be solved exactly to 
find: 
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where we have used the notation: Hoit) = do(t)/a,o(t), 
tt A = A/3H^(t ) and il M = K,p M0 (t )/3H$ (t ) with the 
condition do (to) = 1- On the other hand, by taking the 
trace of @ in this case, ie. f(R) = and dust, we find: 



R (t) - 4A = KpMo(t) 



(11) 



Now we consider again equation @ but in the case where 
we have an arbitrary function f(R) in the action but not 
DE, i.e. we set A = so that T% = - diag(p M , 0, 0, 0). 
Then we have: 

3[l + f'(R)}--l[R+f(R)}-3-Rf"(R) = -k Pm (12) 
a 2 a 

where we have eliminated the subindex in the different 
quantities to avoid any confusion with the previous case. 
At this stage it is clear that above equation (|12f> solu- 
tions will depend on the function f(R): different choices 
for this function will lead to different evolutions of the 
universe for the same initial conditions. However, our 
approach to the problem will be to find a function f(R) 
so that the solution a(t) of (|12fl will be exactly the same 
as the solution in l|10l) that we got by using standard cos- 
mology and which fits the present cosmological data. In 
other words we want to find a f(R) such that: 



a(t) = a (t) 



(13) 



for the same initial (or better present, i.e. t = to, con- 
ditions). If it were possible to find this function f(R) 
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then it would be possible to avoid the necessity for intro- 
ducing any cosmological constant just by modifying the 
gravitational sector of the action. In the following we will 
show that such a function happens to exist and we will 
give its precise form. In order to do that we first notice 
that having a(t) = <io(t) in this period of universe life 
clearly implies R(t) = Ro(t) and then we can substitute 
R by Rq in ljl2|l . On the other hand we will write the 
matter density as the former matter density plus a new 
contribution, ie. pM{t) — PMo(t) + Ap(t). 

Assuming that matter for arbitrary f(R) is still non- 
relativistic (i.e. dust) in this cosmological era we have 



Ap{t) = Ap(t ) 



apfto) 
a (t) 



where according to JfjJ) particularized for a 
(flip we can write 



Qpfto) 
ao(f) 



and then l|14|l becomes 



Ap(t) = -r) 



R a -4A 

KpMo(to) 



Ro-iA 



(14) 
ao, and 

(15) 
(16) 



where we have introduced the parameter r\ = 
— Ap(to)/pMo(to) so that matter density is written as 
PM(t;r]) = (1 — T))pMo(t)- Finally the last term on the 
l.h.s. of H12|) can be written in terms of the scalar cur- 
vature by differentiating (fill) and using (JSJ. Thus we 
get 



3 [R - 3 A] [R - 4 A] f"(R ) 
- lf(R )-A- V [R Q -AA] 



-~Ro + 3A 







f'(Ro) - 

(17) 



This last equation can be considered as a differential 
equation for the function f(R) (in the following we will 
omit the subindex in R since no confusion is possible). 
(|17f) is a second order linear equation so two initial condi- 
tions are needed to solve it: /(0) and /'(0) for instance. 
The natural choice for these initial conditions will be the 
following: Firstly we do not want to have any cosmolog- 
ical constant in our action, so that f(0) — 0. Second- 
ly we want to recover the standard EH action for low 
scalar curvatures without redefine the Newton constant, 
i.e. /'(0) = 0. Moreover we want f(R) to be an ana- 
lytical function at the origin so that R — should be a 
solution for the field equations in vacuum. This is an ex- 
tremely important requirement since it allows Minkowski 
and Schwarzschild to be vacuum solutions. 

With these initial conditions (|17|l can be solved by us- 
ing standard methods. A particular solution is: 



The homogeneous equation associated with (|17fl is a 
Gauss equation solved in terms of hypergeometric func- 
tions. The general solution of the homogeneous equation 
can be written as: 



f h (R) = A[K+f+(R) + K-f-(R)] 



where 



f± (R) = a a± 2 F\ [a± , 1 + a± - c; 1 + a± - a T ; - 



with a = 3 — R/A and 



a± 



(7± V73) 



-1/2 



(19) 



(20) 



(21) 



The 7] dependent constants K + and K_ must be deter- 
mined from the initial conditions given above. Numer- 
ically we find: K+ = 0.6436(-0.9058?y + 0.0596) and 
K- = 0.6436(-0.2423?7+ 3.4465). 

The hypergeometric functions given in l|20|) are gen- 
erally defined in the whole complex plane. However we 
want to have a real gravitational action. In principle 
this is very easy to achieve since the coefficients in the 
equation Ijl7|l and the constants K + and K- are all of 
them real. Then it is obvious that the real part of (|2C)|) 
is a proper solution of homogeneous equation associat- 
ed with i|17|) . Thus the function we are seeking can be 
written as: 



f(R) = f p (R) + Re[f h (R)} 



(22) 



f P (R) 



-rjR + 2A[rj- 1] 



(18) 



Nevertheless situation is a bit more complicated. The 
homogeneous equation has three regular singular points 
at i?i = 3 A, i?2 = 4 A and R3 = 00. This results 
in the solution fh{R) having two branch points R\ and 
i?2. More concretely there are two cuts along the real 
axis: one from minus infinity to R\ and another from 
i?2 to infinity. Thus one must be quite careful when 
interpreting 1)22(1 . From minus infinity to R\ there is 
only one Riemann sheet of fh(R) where /(0) and /'(0) 
vanish and therefore this is the one that we have to use 
to define f(R). From Ri to R2 the real part of fh(R) 
is well defined. Finally from R2 to infinity there is only 
one Riemann sheet producing a smooth behavior of f(R). 
To reach this sheet one must understand R in the above 
equation as R + ie. 

At the present moment we do not know if this ana- 
lytical structure has any fundamental meaning or it is 
just an artefact of our construction. Much more impor- 
tant is the fact that the function R + f p (R) + fh{R), 
which is the analytical extension of our lagrangian, is 
analytical at R — 0, having at this point the local behav- 
ior R + 0(R 2 ). Therefore our generalized gravitational 
lagrangian R + f{R) does guarantee that R = is a 
vacuum solution whilst it reproduces the current evolu- 
tion of the universe without any cosmological constant. 
Once we have obtained the f(R) function it is possible 
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to check out our result by solving l|12|) in terms of a(t) 
for the f(R) given in (|22|l . This is done by rewriting 
equation (|12fl in terms of a(t) by using © and (JSJ) to- 
gether with the pM definition in terms of 77. This pro- 
cess requires three initial conditions given for the present 
time to which are: a(to) — l,a(ta) — Ho{to)a(to) and 
a(io) = —<lo a (to)Ho(to) with go — —0.61. The numerical 
solution a(t) shows a nice agreement with ao(t) given in 
Thus we can be sure that our gravitational action 
proportional to R + f(R) provides the same cosmic evo- 
lution as EH action R — 2A in a dust matter universe. 
Therefore, our model will verify, in the same range of 
precision, all the experimental tests that the standard 
cosmological model verifies in the present era. 

Notice also that in principle this can achieved for any 
value of 77, i.e. for any desired amount of matter. Nev- 
ertheless some restrictions should be imposed over the 
parameter 77. For instance it is obvious that in a dust 
matter dominated universe pM(t',v) ^ implies 77 < 1. 

Much more stringent bounds can be set on 77 by de- 
manding our model to work properly back in time up to 
Big Bang Nucleosynthesis (BBN) era. Observations indi- 
cate that the cosmological standard model fits correctly 
primordial light elements abundances during BBN with a 
10% of relative error for Ho (t) . Therefore by the time of 
BBN, departure of our model from the standard cosmol- 
ogy must not be too large and <(T2")l should give similar 
behavior to the one given by the standard Friedmann 
equation © where now p M0 = Po ust (t) + p r Q adlatlon (t) . 

At BBN era cosmological constant is negligible com- 
pared with dust and radiation densities. The scalar cur- 
vature is of order 10 -39 eV 2 (with h = c = 1 for these 
calculations) and by that time dust and radiation den- 
sities are of the order of 10 16 eV A and 10 21 eV 4 respec- 
tively. Since R ~ i?o we can rewrite Ijl2|l as a modified 
Friedmann equation as follows 

H 2 (f] = R 2 (f] [ 10 5 R-yR+^Rf'(R)-f(R)] ' 

U oU [l0 5 R[l + r{R)-3f"(R)(l- V )R] ' 

(23) 

As it was commented above, to reproduce light elements 
abundances it is required that H 2 (t) = if, 2 (/j)(l±0.2) for 
curvatures of order Rb b n ■ This implies that the second 
factor on the r.h.s. of 11' Ml should be between 0.8 and 1.2 
by that period. Thus in order to match our f(R) gravity 
model with the standard cosmology at the BBN times 
we need to tune 77 to a value about 0.065. Therefore the 
matter content of our model is not too different from the 
one in the standard cosmology and the difference is in 
fact smaller than experimental precision in 0. 

To conclude we have succeeded in finding a f(R) grav- 
ity which exactly reproduces the same evolution of the 
universe, from BNN up to the present time, as stan- 
dard cosmology, but without the introduction of any 
form of DE or cosmological constant. The gravitation- 
al lagrangian R + f(R) is analytical at the origin and 



consequently R = is a vacuum solution of the field 
equations. Therefore Minkowski, Schwarzschild and oth- 
er important R = GR solutions, with A = 0, are also 
solutions for this f(R) gravity. The price we have to 
pay for all those good properties is that our lagrangian, 
considered as a function of R, has a very complicated 
analytical structure with cuts along the real axis from 
infinity to R = 3A and from R = 4A to infinity. Ob- 
viously the only reasonable interpretation of our action 
is as some kind of effective action. In classical physics 
one typically starts from some action principle, obtains 
the corresponding field equations and finally solves them 
for some initial or boundary conditions. In this work we 
have proceeded in the opposite way: we started from so- 
lutions obtained in the standard cosmological model and 
then we have searched for an action that, possessing cer- 
tain properties, gives rise to field equations having the 
same solutions. Classical actions are of course real but 
effective quantum actions usually have a complex struc- 
ture coming from loops and related to unitarity. The 
presence of an imaginary part in the action, evaluated 
on some classical configuration, indicates quantum lost 
of stability by particle emission of this configuration 
Therefore it is tempting to think that our action could 
have some interpretation in terms of an effective quan- 
tum action. However, our action determination proce- 
dure does not allow to make such a kind of statement. 
Complications in the action could be just an artifact of 
our construction. In any case we have shown that such 
action exists and it reproduces the present universe evo- 
lution without DE having R = as a vacuum solution. 
The complicated structure of this action may be an in- 
dication that the cosmological constant problem is even 
much harder to solve than we have previously thought. 
Obviously much more insight and research are needed in 
order to get further progress in this issue. 

Acknowledgements: We would like to thank J.A.R. 
Cembranos and very specially A.L. Maroto for useful 
comments and discussions. This work has been partial- 
ly supported by the DGICYT (Spain) under the project 
FPA2005-02327. 



[1] A. G. Riess et al. [Supernova Search Team Collaboration], 

Astrom. J. 116, 1009 (1998) arXiv:astro-ph/9805201 . 
[2] S. Perlmutter et al. [Supernova Cosmology Project 

Collaboration] , Astrophys . J. 517, 565 (1999) 

|arXiv:astro-ph /98 12133] . 
[3] J. L. Tonry et al. "A strophys. J. 594, 1 (2003) 

[arXiv:astro-ph/0305008l . 
[4] D.N.Spergel et al. Wilkinson Microwave Anisotropy 

Probe (WMAP) Thre e Year Results: Implications for 

Cosmology, astro-ph/0603449 2006. 
[5] S. M. Carroll, A. De Felice, V. Duvvuri, D. A. Easson, 

M. Trodden and M. S. Turner, Phys. Rev. D 71 (2005) 

063513 arXiv:astro-ph/0410031 . 



5 



[6] A. A. Starobinsky, Phys. Lett. 91B, 99 (1980). 
[7] J. D. Barrow and A. C. Ottewill, J. Phys. A 16, 2757 
(1983). 

[8] A. Dobado and A. L. Maroto Phys. Rev. D 52, 1895 
(1995). 

[9] S. M. Carroll, V. Duvvuri, M. Trodden and 
M. S. Turner, Phys Rev. D 70 (20.04) 043528 
arXiv:astro-ph/0306438 . 



[10] G. J. Olmo, Phys. Rev. Lett. 95 (2005) 261102 

arXiv:gr-qc/0505101 . 
[11] J. A. R. Ce mbranos, Phys. Rev. D 73 (2006) 064029 

arXiv:gr-qc/0507039 . 

[12] T. P. Sotiriou and S. Liberati arXiv:gr-qc/0604006]. 
[13] A. Dobado and A. L. Maroto, Phys.Rev.D60, (1999) 

104045 



